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Abstract. A coagulation process is studied in a set of random masses, in which two 
randomly chosen masses and the smallest mass of the set multiplied by some fixed 
parameter uo € [—1,1] are iteratively added. Besides masses (or primary variables), 
secondary variables are also considered that are correlated with primary variables 
and coagulate according to the above rule with uj = Q. This process interpolates 
between known statistical physical models: The case a; = — 1 corresponds to the 
strong disorder renormalisation group transformation of certain disordered quantum 
spin chains whereas to = 1 describes coarsening in the one-dimensional Glauber-Ising 
model. The case w = is related to the renormalisation group transformation of a 
recently introduced graph with a fat-tail edge-length distribution. In the intermediate 
range — 1 < < 1, the exponents and (3^ that characterise the growth of the 
primary and secondary variable, respectively, are accurately estimated by analysing the 
differential equations describing the process in the continuum formulation. According 
to the results, the exponent varies monotonically with lo while j3uj has a maximum 
at w = 0. 
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1. Introduction 

Coagulation processes arise in various areas of physics; one may think of polymerisation, 
growth of ordered domains in non-equilibrium magnetic systems p], dynamics of 
droplets when water condenses on non- wetting surfaces [2], etc. The substance, 
or "mass" that aggregates is very frequently not conserved during the process: for 
example, agglomerating insoluble inclusions in molten metal may be lost from the melt 
by attachment to the wall of the vessel [3J. Therefore the theoretical investigation 
of the kinetics of such non-conserving coagulation processes is of great importance. 
Moreover, the models developed for the description of such systems may show interesting 
behaviour: the Smoluchowsky equation with certain coagulation kernels exhibits 
gelation transition and, in general, even the simplest models with conserved mass 
may have non-trivial solutions, see e.g. Ref. [3] and references therein. Beside quite 
realistic ones, there is a special class of (possibly non-conserving) coagulation models 
where only the actually smallest one among the masses is active while the other masses 
are temporarily inert. This type of extremal dynamics can be regarded as a rough 
approximation for models where the reaction rates are decreasing functions of the mass 
of particles. In what follows, we shall survey three processes with extremal dynamics 
in detail. We mention, however, that models of this type have also been introduced 
in the context of dynamics of growing and coalescing droplets |2] or multispecies pair 
annihilation reactions [1]. 

In the one-dimensional Glauber-Ising model started from a random initial state at 
zero temperature, the domain walls move as independent random walkers and annihilate 
upon meeting. While the closest pairs of walls come together and annihilate, the 
other domain walls hardly move. A simplified model of evolution of distances Xj 
between adjacent walls can be formulated as follows [Sj [6]. The shortest interval 
Xm is eliminated together with the two adjacent intervals Xi and X2 and replaced 
by X = Xi + X2 -|- Xm- As the density of walls tends to zero, the distributions of 
intervals at different times become self-similar, depending on a single time-dependent 
length scale, and the corresponding scaling function can be calculated exactly [5l[7l[8]. 
Another quantity of interest is the fraction of space which has never been traversed by a 
domain wall. The length Yi of such parts of intervals transforms in the way Y = Y1 + Y2 
when the shortest interval is eliminated. The characteristic value of X depends on 
the fraction c of the initial intervals that have not yet been eliminated as X ~ c~", 
obviously, with a = 1, while it has been found that Y ~ c~^, where the persistence 
exponent (3 = 0.82492412 ... is the zero of a transcendental equation [7]. In addition to 
this, the autocorrelation exponent has also been exactly calculated in this model [9J. To 
obtain this quantity, the overlap Zj of an interval with its initial state that transforms 
as Z = Zi + Z2 — Zm had to be considered. Later, a generahsation of persistence 
has been studied in the same model, which required to introduce an auxiliary variable 
transforming a.s Y = Yi + Y2 + pY^ [lOj. Here, the generalised persistence exponent 
has been found to vary monotonically with the partial survival factor p in the range 
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-1 < p < 1. 

The next example is the strong disorder renormahsation group transformation 
of inhomogeneous quantum spin chains [H]. Here, the degrees of freedom related to 
the largest coupling (a bond between neighbouring spins or a local external field) are 
eliminated one after the other. In terms of logarithmic couplings, Xi, the renormahsation 
rule generally reads as X = Xi + X2 — X^, where X is a newly formed effective variable 
and Xi,X2 are variables adjacent to the smallest one, X^- For the relation between 
these variables and the couplings in the particular Hamiltonians we refer the reader to 
Ref . |12j . A variable Yi that transforms according to the rule Y = Y1 + Y2 under such a 
renormahsation step can be interpreted in the case of a particular model, the transverse 
field Ising chain, as the magnetic moment of a spin. For this process with i.i.d. random 
initial variables Xi, which corresponds to critical spin chains, the distribution of X flows 
again to a fixed point where it shows scaling behaviour. The characteristic value of X 
increases in the course of the process as X ~ c~" with a = 1/2, while the variable Y 
grows as y ~ c"^ with (3 = {1 + ^/5)/A = 0.809016 . . . [I2]. Note that the coagulation 
rules in the above two models differ only in the sign of X^, which leads to different 
exponents a and p. 

Our third example is a random graph where three edges emanate from each node, 
and which is built on a regular one- dimensional lattice by adding long edges in the 
following way. To each edge of the one-dimensional lattice that we call short edges, a 
random weight Xj is assigned. Defining the length of a path as the sum of weights of 
the edges it contains, the closest pair of nodes of degree 2 with respect to this metric is 
chosen and connected by an edge of unit weight. This step is then iterated until all nodes 
become of degree 3 [13]. For this graph, a renormahsation procedure can be formulated 
where loops are eliminated step by step in reversed order compared to the construction 
procedure. Formally, the short edge with the minimal weight X^ is eliminated together 
with the nodes it connects, as well as with the neighbouring short edges with weights 
Xi, X2 and a new effective short edge is formed with a weight calculated asymptotically 
as X = Xi + X2. According to numerical results, the characteristic value of effective 
weights grows as X ~ c~" with a = 0.826(1) |13]. This exponent characterises at the 
same time the diameter of finite graphs with N nodes with respect to the above metric 
via D{N) ~ X". 

As can be seen, these seemingly different problems can be treated in a common 
framework and can be interpreted as coagulation processes with extremal dynamics. In 
the first example, the total sum of the variables Xj is conserved while in the latter two 
cases it is not. We will study in this work a coagulation model controlled by a parameter 
u that interpolates continuously between the first two models and incorporates the 
third one as a special well. We are interested in the exponents a^^ and fS^j for 

intermediate values of the parameter u and shall provide accurate estimates for that 
is obtained as the root of a transcendental equation while is accurately determined 
by the numerical analysis of a system of non-linear differential equations. We shall 
see that varies monotonically between the corresponding values of the two marginal 
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models, while, unlike the generalised persistence exponent of the model with partial 
survival mentioned above [TO] , the exponent jS^^ shows a maximum when uj is varied. As 
can be seen, the transformation rule of the variable Y does not depend directly on the 
parameter lv but it is influenced indirectly via the correlations emerging between X and 
Y, the strength of which is controlled by lj. Therefore our results may contribute to the 
understanding of the role of correlations in such models. Moreover, these investigations 
provide an accurate estimate for the diameter exponent of the graph quoted above, for 
which we obtain a = 0.82617561 in agreement with the previous numerical result. 

The rest of the paper is organised as follows. In Section [21 the model and its 
continuum description is introduced. In Sections [3] and HI the way of approximative 
determination of the exponents and [3^ is presented. Some calculations are given in 
the Appendix. Finally, results are discussed in Section [5l 

2. The model and its continuum formulation 

2.1. Definition of the model 

Let us consider a finite set of positive vectors Vi = (Xj, Yi) indexed by the 
integers i = 1,2, . . . , N. We assume, moreover, that is odd. The vectors 
are independent, identically distributed random variables drawn from a continuous 
distribution p{X, Y)dXdY, for which we require that all moments exist. The first 
components Xj and the second components Yi are called primary and secondary 
variables, respectively. Assume, furthermore, that uj G [—1, 1] is a fixed real number. 
Now, the following procedure is considered on this set. The vector Vm with the smallest 
primary variable is chosen and, at the same time, two further vectors Vi and Vj are 
chosen at random from the set. These three vectors are removed and a new vector V 
with components 

X = Xi -\- Xj -\- LuXjji 

Y = Y, + Y, (1) 

is added to the set. Thereby the number of vectors in the set is reduced by two. Note 
that the vectors remain independent after such an operation and that 

X > Xi, Xj, Xjn (2) 

even for uj = —1. This step is then iterated until a single vector Vn = {Xn,Yn) is 
left in the set. In this general formulation, the cases u = 1,-1,0 correspond to the 
three models in the order as they were quoted in the Introduction. Based on the known 
asymptotical behaviour of Xn and Yn for large N in the marginal cases u = —1, 1, we 
expect 

Xjv-AT"- and Y^v ~ (3) 

to hold also for intermediate parameter values — 1 < u; < 1 with some exponents and 
that may depend on u. 
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2.2. Continuum formulation 

Now, we consider the continuum limit N oo and introduce the probabihty density 
Pr{X) of the primary variable that has the support F < X < oo and that depends 
on the lower boundary F as a parameter. The function PriX) is normalised as 
Pr{X)dX = 1 for any F. Following Ref. [7j, we consider, furthermore, the expected 
value Yr{X) of the secondary variable under the condition that the primary variable is 
X. In the continuum limit, the system is described by these two functions of X, which 
depend on the lower boundary of the support F as a parameter. The inequality ([2]) 
implies that, as the fraction of vectors cr that have not yet been eliminated decreases in 
the course of the coagulation process, the lower edge F of the distribution continuously 
increases. As it is shown in the Appendix, one may write the following differential 
equation for Pr{X): 

-^r^ = Pr(r)e[X-(2+.;)F] Pr{X')Pr{X-X' -ujT)dX\ (4) 

where Q{X) is the Heaviside step function. The fraction cr is related to F as 
dcr/cr = — 2Pr(r)(iF or, equivalently, 

^ - -2PHT)cr. (5) 

The function Qr{X) defined as 

Qr{X) ^ PriX)YriX), (6) 
can be shown to obey the differential equation 



2Pr(r)e[X-(2+cj)F]^ Qr{X')Pr{X-X'-ujr)dX'.{7) 



dT 

The derivation of this equation is given again in the Appendix. 

2. 3. Fixed point solution 

In the marginal cases uj = —1, 1, it is known that, for any well-behaving initial 
distributions p{X, Y) with finite moments, the solutions of Eqs. (jlj) and (JTj) tend to 
a universal fixed point solution P^{X), Q^{Y) in the limit F ^ cxo that has the scaling 
property 

P*{X) = F-V(X/F) 

QUX)=T'-''giX/T), (8) 
with some number 5^ that is related to the growth exponents a^ 

Sc, = Poj/Oiuj- (9) 

Therefore we expect this to hold also for intermediate parameter values — 1 < < 1 
with some (a priori unknown) exponent 6^^ that may depend on u. Indeed, the functions 

t This can be seen from the equation Fr(r) = Q^{T)/P^{T) = T^'- g{l)/ f{l) that indicates the 
asymptotical relation Y ~ X^" between the typical values of primary and secondary variables. 
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in Eq. solve Eqs. (jl]) and ([7]) provided that the universal scaling functions f{x) and 
g{x) satisfy the following differential equations: 

d[xfix)] ^ _f^Q^^ - 2 - cu) r^'^" fix') fix -x'- u)dx' (10) 



dx 

d[x^-^-gix)] rx-l-u. 



-X 



/X— L—LJ 
gix')fix-x' -iu)dx', (11) 

where the notation /i = /(I) has been used. For an alternative derivation of these 

equations in the case u = 1, see Ref. [7j. Using the fixed point solution, Eq. can be 

integrated yielding the asymptotic relation in the large F limit: 
]_ 

F-Cr'-'^ (12) 
Comparing this with Eq. ([3]), we obtain the relation: 

(^00 = 7^- (13) 
3. Approximative determination of 

As can be seen, Eq. (fTOl) does not contain gix) and together with Eq. (fT3l) it constitutes 
an autonomous problem for the calculation of the exponent a^j- For the special case 
u = —1, the solution of Eq. ffTOl) is of simple form: fix) = e~^~^^; this yields a_i = |. 
In the other marginal case, u = 1, the Laplace transform of the solution is known 
E] and tti = 1. In the case — 1 < u; < 1, where Eq. flTU]) is not soluble, we shall 
construct an approximative solution that enables us to give an accurate estimate of 
a^. An alternative way related to the numerical analysis of the Laplace transforms is 
presented in the next section. 

Some properties of the scaling function fix) can be easily established by 
investigating Eq. ( ITOl) without knowing the exact solution. Apparently, the r.h.s. of Eq. 
(fTOj) and, as a consequence, fix) is non-analytical at x = Xi = 2 + u. But, as fix) itself 
appears on the r.h.s. as a convolution with a shifted argument x — 1 — u, the r.h.s. as 
well as fix) must be non-analytical also at a; = X2 = Xi + l + u. Iterating this argument, 
it turns out that there are infinitely many points where fix) is non-analytical. To be 
precise, one can show by recursion that the 2nth derivative of fix) is discontinuous a.^ 

Xn = I + il + uj)n, n = 0,1,2,.... (14) 

Furthermore, the function value of fix) at some x' is determined by fix) in the restricted 
domain il,x' — 1 — u). Due to this property, fix) can be constructed in the intervals 
[x„,x„+i] step by step starting with n = 0. However, the solution is more and more 
complicated for increasing n as it contains multiple integrals that cannot be evaluated 

§ For a more direct way to this result in the case w — 1, where the cxphcit form of the scahng function 
f{x) is available, see Ref. [S]. 
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analytically. In the domains x^+i], n = 1, 2, . . ., the function f{x) can be written in 
the following form: 

1 " 

fix) = - E /f ^'C'f +i)(a;), xn<x< x^+u (15) 

X 1=0 

whereas f{x) = if x < Xq- The functions Cl^^~^^\x) are independent of /i and the first 
three of them read as 

CW(x) = 1, (16) 

Substituting an exponential trial function in Eq. (fTOi) . we obtain that the asymptotical 
solution foo{x) in the limit a; — > cxd is of the form 

/oo(x) = ^e-<^^-\ (17) 
Ji 

with the some number a that is not fixed by this substitution. The graph of f{x) for 
d; = is shown in FigUl As can be seen, /(x) tends rapidly to fooi^) for increasing 




Figure 1. Graphs of the functions f{x), g{x) and foo{x) for = 0. The former two 
are obtained by numerical integration of Eqs. PT!]) and ITTI) . whereas the latter is given 
in Eq. dni). 



X. This suggests an approximation for f{x) in which f{x) is replaced by the simple 
asymptotical function foo{x) for large x. To be precise, the nth (n = 0,1,2,...) 
approximant /'-"^(x) is defined as 

/(")(x) = /(x), if x<x„ 

/W(a;) = /^(x), otherwise. (18) 
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The two unknown parameters /i and a are determined by the requirements that /'■"^(x) 
is continuous at x = x„, i.e. 



and that it is normahsed as 

f{x)dx + / foo{x)dx 



Xo 



(19) 



(20) 



Using the expression in Eq. (ITS!) . straightforward calculations result in that the nth 
approximant /i"^ {n > 0) is the root of the following transcendental equation: 



i=0 



+ 



Xn + OJ 



n-1 



j=0 



n-l 



i=0 



where the function A^^^*"'"^-' {x) has been introduced as 

We have numerically calculated the root of Eq. fl?Il) and the nth approximant a^"-* of 
Oo; by using Eq. (fT3l) for n = 1,2,3 and for several values of uo. This has necessitated 
the numerical evaluation of the integrals in Eq. (!22|) for n > 1. Results are shown in 
Fig. [2] and some numerical values are given in Table I. As can be seen, the approximants 
Q,£"') converge rapidly with increasing n and they increase monotonically with oj. The 
best estimate for the diameter exponent of the graph cited in the Introduction is 
= 0.82617561. 
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Figure 2. The third approximant al, of the exponent plotted against u. 
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0c = 5^a^^^ 


-0.9 


0.54752760 


0.54752815 


0.54752815 


1.48973578 


0.81567227 


-0.8 


0.59036797 


0.59037862 


0.59037860 


1.38841226 


0.81968889 


-0.7 


0.62906729 


0.62911723 


0.62911708 


1.30687751 


0.82217896 


-0.6 


0.66421085 


0.66434418 


0.66434376 


1.23995279 


0.82375490 


-0.5 


0.69632781 


0.69659263 


0.69659189 


1.18399594 


0.82476197 


-0.4 


0.72586754 


0.72630756 


0.72630667 


1.13643762 


0.82540221 


-0.3 


0.75320237 


0.75385266 


0.75385199 


1.09543864 


0.82579860 


-0.2 


0.77863838 


0.77952410 


0.77952421 


1.05965756 


0.82602872 


-0.1 


0.80242716 


0.80356405 


0.80356560 


1.02809664 


0.82614309 


0.0 


0.82477635 


0.82617193 


0.82617561 


0.99999999 


0.82617561 


0.1 


0.84585830 


0.84751339 


0.84751989 


0.97478484 


0.82614953 


0.2 


0.86581708 


0.86772721 


0.86773715 


0.95199472 


0.82608118 


0.3 


0.88477397 


0.88693068 


0.88694457 


0.93126697 


0.82598219 


0.4 


0.90283179 


0.90522369 


0.90524198 


0.91230963 


0.82586098 


0.5 


0.92007834 


0.92269199 


0.92271501 


0.89488491 


0.82572374 


0.6 


0.93658910 


0.93940970 


0.93943769 


0.87879701 


0.82557503 


0.7 


0.95242937 


0.95544130 


0.95547441 


0.86388324 


0.82541834 


0.8 


0.96765597 


0.97084323 


0.97088154 


0.85000684 


0.82525595 


0.9 


0.98231868 


0.98566518 


0.98570869 


0.83705175 


0.82508918 


1.0 


0.99646128 


0.99995110 


0.99999976 


0.82492447 


0.82492427 



Table 1. Approximants of the exponents a^^, and for different values ofuj. 



4. Numerical determination of (3^^ 

Next, we turn to the determination of the exponent (and, at the same time, (3^^ 
through Eq. (I9l)), which requires the analysis of the full problem, i.e. the system of 
differential equations (fTOj) and (ITT]) . Prior to this, a few remarks concerning the scaling 
function g{x) are in order. First, as a consequence of the definition in Eq. ([6]), g{x) 
apparently inherits the singularity properties of f{x) discussed in the previous section. 
Furthermore, it can be written in a form analogous to Eq. f|T5l) . In the domain [xq, Xi], 
it has a simple form: 

g{x) = g{l)x^"-\ Xo<x < xi. (23) 

Second, the differential equation (ITT]) gives the scaling function g{x) only up to a 
multiplicative constant. This non- universal constant depends on the initial distribution 
p{X, Y) and it is fixed in a non-trivial way by the original equations (jlj) and ([7j) that 
are valid for any F. Third, the equation ( ITTi) contains the a priori unknown parameter 
6^ that must be fixed by physical considerations about the solution that depends on 6^. 
Namely, the physically acceptable solution must be nonnegative and must have the only 



A non- conserving coagulation model with extremal dynamics 10 

reasonable asymptotics allowed by Eq. (ITTi) : 

goo{x) ~ const ■ xe"""", (24) 

where the number a is the same as that appears in Eq. ( |T7|) . Numerical analysis of Eq. 
( ITTj) shows that these requirements are fulfilled only for a single value of the parameter 

Following Ref. [7], it is, however, simpler to analyse the Laplace transform of the 
equations ffTOl) and ffTTl) . Introducing the functions 



= e-P''f{x)dx, ilj{p) = e-P''g{x)dx, (25) 
the equations (ITOl) and (ITTl) transform to 

P0'(p) = /i[e-">'(p)-e-^], (26) 
p^'(p) = -5^ij{p) - gie-^ + 2f,e-'^P^{p)^{p), (27) 

where the prime denotes derivation by p and gi = g{l). These equations are not soluble 
in the parameter range — 1 < < 1 but asymptotical expressions of the solution can be 
established. The functions (j){p) and ip{p) have the small-p expansions: 

oo oo 

Hp) = E ^(P) =9iT. M"- (28) 

n=0 ra=0 

Substituting these into Eqs. (l26l) and (l27j) . we obtain that the expansion coefficients for 
-1 < uj <1 are given by = 1, 6o = 2/^-^ ^^'^ t)y the following recursion relations 



for n > dJi]: 



^ r). + X]o<i,i,A;<n; i+j+k=n ^ ji ^ 

f -2 

( 



(29) 



i! V2/i-<5cj 2fi ) ^ 2fi-5uj ^ ^0<i,j,k<n; i+j+k=n i\ 



(30) 

±1 — I 

2/1 

Next, we discuss the large-p behaviour of ip{p). The differential equation fl271) can 
have two kinds of asymptotical solutions depending on the parameter 5^. If the second 
term on the r.h.s. dominates, we obtain 

^'{p) ^ -91 — , (31) 
P 

while, if the first term dominates, we obtain 

iIj{p) ~ const ■ p~^'^. (32) 

On the other hand, it follows from Eq. (!25|) that ipip) must have the large-p asymptotics: 
ip{p) ~ 5'ie~^/p[l + 0{l/p)]. Using that the function g{x) is explicitely known in the 
domain 1 < x < 2 + uj (see Eq. (I23l) ). we can obtain a more accurate asymptotical form 

II These series expansions are also valid for uj = —1 with a2 — 5/2, and for lo ~ 1 with ai — —2e^ [7], 
where 7 is Euler's constant, given by 7 = - Inte-^dt = 0.577215 . . .. 



A non- conserving coagulation model with extremal dynamics 



11 



for 1/p ^ 1+uj. Replacing g{x) by the function in Eq. (125]) in the entire domain x > 1, 
the integral in Eq. fl25|) can be evaluated, yielding 



VUp)-*.e-^g(*";^)^, (33) 

This shows that the physically acceptable asymptotics is that given in Eq. fl3T]) whereas 
that in Eq. ( l32l) is non-physical. Numerical analysis of the differential equations ( l26l) 
and (1271) with the correct value of /i shows the following behaviour of the solution when 
the parameter S^j is varied: For small enough S^^, the function ip{p) is non-monotonous 
and tends to zero from below for increasing p as given in Eq. (l32l) : for large enough 6uj, 
the function iplp) decays monotonically to zero again with the asymptotics given in Eq. 
( 132|) . These parameter regimes are separated by a "critical" value of 6i^. At this value, 
the solution decays monotonically to zero with the physically acceptable asymptotics 
given in Eq. (l3Tl) . 

The numerical estimation of 5^^ is based on this scenario: The differential equations 
(126|) and (127|) are integrated from p = to some large p and the true value of is 
selected by the condition that ipi^p) has the correct asymptotics. We have assumed here 
that /i is already at our disposal. This can be obtained either by the approximative 
procedure described in the previous section or, analogous to the above method, from 
the condition that (j){p) has the correct asymptotics given by Eq. (l33l) with 5^ = 0. 

Before presenting numerical results on 5^, we show that, in the case u = 0, the 
assumption on the uniqueness of the value 6uj that corresponds to the correct asymptotics 
implies that 6^ = 1. For u = 0, the primary and the secondary variables coagulate 
according to the same rules, see Eq. ([T]). If these variables are initially perfectly 
correlated, i.e. Xi = hYi with some common constant h for all z, it is obvious that 
«o = /?o- Nevertheless, this equality holds for general initial distributions, as well. 
Indeed, it is easy to check that for = 0, the function 

m = -T^'ip) (34) 

solves Eq. ( 1271) provided that 5q = 1 and 0(p) is the solution of Eq. ( l26l) . In terms of 
the scaling functions, equation (1M|) reads as g{x) = xf{x)^. As the asymptotics of the 
solution in Eq. (l34l) is physically acceptable, we conclude that 



So = 1. (35) 
The details of the numerical determination of are the followings. The best 

/o\ _____ 

approximant that we have, fl , has been substituted in Eqs. ( l26l) and ( 1271) and a trial 
value for 6^ has been chosen. As the derivatives 0'(p) and ip'{p) calculated from these 
equations are of the form 0/0 at p = 0, the functions (f){p) and ipi^p) were first calculated 
at p = 0.05 by using the small-p expansion in Eq. (l28l) in order to avoid numerical 
uncertainties of the integration in the vicinity of p = 0. Then, starting from p = 0.05, 
the differential equations were integrated by the Bulirsch-Stoer method [H] to some pf, 
where ipi^Pf) is compared to the asymptotical form in Eq. (l33l) . In practice, we have 
monitored the derivative ip'{pf) rather than ip{pf) and pf = 12 was sufficiently large so 
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that the asymptotical value (at the critical 6^^) is reached within the numerical accuracy 
of the integration. The true 6^ was then selected by the condition ip'ipf) = i^'ooiPf)- The 
exponent 5^^ determined in this way is plotted against cu in Fig. [3] whereas P^i = S^ja^J^^ 
is plotted against u in Fig. |H Some numerical values can be found in Table [TJ As can 
be seen, 6^^ decreases monotonically with u but (3^0 has a maximum at u; = 0. The latter 
observation will be explained in the next section. 
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Figure 3. Numerically calculated exponent plotted against uj. 
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Figure 4. Numerically calculated exponent Puj — <5cjaL plotted against lu. The 
horizontal lines indicate known values in the marginal cases w = — 1, 1. 

5. Discussion 

We have shown that two problems, the renormalisation group procedure of certain 
disordered quantum spin chains and a simple model describing coarsening in the 
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Glauber-Ising model can be deformed into each other by varying a single parameter. The 
interpolating model is a special type of non-conserving coagulation process where only 
the actually smallest mass is active. In the range — 1 < < 1, the exponents a^j and P^j 
that characterise the growth of the primary and the secondary variables, respectively, 
vary continuously with u. The latter exponent exhibits a maximum, which contrasts 
with the monotonous dependence of the generalised persistence exponent on the partial 
survival factor that appears directly in the transformation rule of the secondary variable 
|10] . Although, we have focused on the range — 1 < < 1, the equations written down 
in this work are valid also for u > 1. In that case, the growth of the primary variable 
becomes super- linear, meaning that > 1. 

An intriguing feature of the process studied in this work is the universality with 
respect to the initial distribution of the variables: For a fixed u, any sufficiently rapidly 
decaying initial distribution tends at late times to a universal distribution that displays 
scaling. Although, the process is universal in this sense, we have pointed out that it is 
sensitive to the variations of the reaction rules parameterised by uj. The dependence 
of on UJ is obvious since the transformation rule of the primary variable contains u 
explicitely. The growth of the secondary variable is, however, affected by in a more 
subtle way. Focusing on the secondary variables, the difference to the process of primary 
variables with u; = is that, here, not exactly the smallest variable is removed from the 
set. This is the reason for that is unequal to ao for uj ^ 0. Nevertheless, for any 
u, the removed secondary variable is typically relatively small since Xi and Yi become 
positively correlated in the course of the process. Due to these correlations, the strength 
of which is controlled by uj, the variation of (3^^ is relatively slight. Indeed, it is by an 
order of magnitude smaller than that of a^. 

For u = 0, we have shown that ao = even if the primary and the secondary 
variables are initially not perfectly correlated. This can be understood also on a 
microscopic level since, in this case, the vectors in the set are sums of an increasing 
number of initial vectors. Thus, the ratios Xi/Yi tend stochastically to a common 
constant in the limit F — > cxo for all i. In words, the two types of variables become 
asymptotically perfectly correlated for u = 0. Now, we are in a position to understand 
why the exponent f3i^ is maximal at u = 0. At that point, the correlations are (at 
least asymptotically) perfect and almost always the smallest one among the secondary 
variables is removed. For u ^ 0, however, the correlations are no longer perfect and, as 
a consequence, not strictly the smallest secondary variables are eliminated. Therefore 
the fastest growth of Y is realized at u; = 0. 

In a general aspect, the benefit of the analysis carried out in this work is that 
the numerical technique developed here for obtaining accurate estimates of the growth 
exponents may also apply to other non-soluble coagulation processes with extremal 
dynamics. 
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Appendix A. 

When the primary variables in the infinitesimal interval (r,r + dT) are eliminated in 
the course of the process, we may write for the change of the probability density Pr{X): 

Pr+dr{X) = |Pr(X) + Pr(r)rfr J dX^ J dX2Pr{Xi)Pr{X2) 

[6iX - Xi - - ujT) - 5iX - Xi) - 6iX - X^)]} (A-l) 

Here, the first term of the integrand is related to the newly generated primary variable 
while the other two terms are related to the removed ones. The factor — 2Pr(r)o^r] 
ensures that the distribution remains normalised. Expanding the l.h.s. of Eq. ( lA.ip 
and retaining only terms of order dT, we arrive at the differential equation (jlj). 

The expected value Yr{X) of the newly generated secondary variable under the 
condition that the generated primary variable is X is given as 

Yr{X) = (A 
II dXi J dX2[Yr{Xi) + Yr{X2)]Pr{Xi)Pr{X2)6{X - Xi - X2 - ioT)^ /UX) 

2 1 dX'PriX')YriX')PriX - X' - ujT)/UX), 

for X > (2 + iu)T, where the function /^(X) = dX' PriX')Pr{X - X' - uT) 

has been introduced. The expected value Yr+driX) can then be written as the weighted 
average of Y^lX) and Yr{X) as follows: 

Yr+dr{X) = _ (A.3) 

[PrjX) - 2dTPr{T)Pr{X)]Yr{X) + dTPr{T)UX)Yr{X) 
Pr{X) + dVPriV) [UX) - 2Pr{X)] 
This leads to the differential equation 

dPriX)Yr{X) ^ p^(r)e[X - (2 + u;)T]UX)Yr{X), (A.4) 

where we have made use of Eq. (jH). Rewriting this equation in terms of Qr{X) given 
in Eq. ([6]), one arrives at Eq. ([7j). 
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